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Abstract 

We address the possibility of spiral ferrimagnetic phases in the mean-field phase 
diagram of the two-dimensional (2D) Hubbard model. For intermediate values 
of the interaction U and doping n, a spiral ferrimagnetic phase is 

the most stable phase in the {n,U) phase diagram. Higher values of U lead to a 
non-spiral ferrimagnetic phase. If phase separation is allowed and the chemical 
potential /i replaces the doping n as the independent variable, the (/i, U) phase 
diagram displays, in a considerable region, a spiral (for 6 < U/t < 11) and non¬ 
spiral (for higher values of U) ferrimagnetic phase with fixed particle density, 
n = 0.5, reflecting the opening of an energy gap in the mean-field quasi-particle 
bands. 

Keywords: A. Magnetically ordered materials, G. Grystal structure and 
symmetry, D. Electron-electron interactions, D. Electronic band structure, D. 
Phase transitions 


1. Introduction 

The 2D Hubbard model remains the most important open theoretical prob¬ 
lem in the field of the strongly correlated electronic systems, despite all efforts 
fuelled by the advent of the high-Tc superconductivity. [TJ [2] At half-filling, the 
spin dynamics of the 2D Hubbard model is described by the Heisenberg anti¬ 
ferromagnetic exchange term. [3] Away from half-filling, the movement of holes 
through the spin background generates additional spin mixing. The competition 
between the Heisenberg exchange and the spin configuration mixing generated 
by hole hopping in the 2D Hubbard model is still far from understood. m mil] In 
particular, there is no consensus regarding the ground state magnetic phase dia¬ 
gram of the 2D Hubbard model and different authors obtain different mean-field 
(ME) phase diagrams depending on the magnetic phases allowed. [7] Tradition¬ 
ally, one considered ferromagnetism, antiferromagnetism and paramagnetism 
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Figure 1: The 2D lattice and its four sublattices A, B, C and D. We consider two situations: 
(i) rriA = 'rriD = and uib = ttic = (ii) mA = and ttib = tub = 'm 2 - 


phases. laiaiiniiiiiin] The complexity of the MF phase diagram was increased 
with the introduction of spiral phases [13], which appear between the ’’usual” 
magnetic phases in the diagram. This complexity was further increased by the 
consideration of spatial phase separation. nsnsiiini 

In this paper, we extend the results above mentioned, by introducing the 
possibility of a spiral ferrimagnetic phase, that is, a ferrimagnetic phase such 
that the orientation of magnetic moments changes along the lattice (see Fig. 
[^. More precisely, we study the 2D Hubbard model using the Hartree-Fock 
approximation in a square lattice decomposing the lattice in four square sub¬ 
lattices (A, B, C and D as in Fig. and allowing different amplitudes for 
magnetizations of the spiral phases in the sublattices. Note that, even under 
the MF approximation, when four sublattices are considered, it is not possible 
to obtain the analytical form of the spectra of the 2D Hubbard model. Our MF 
approach to the 2D Hubbard model follows that of Dzierzawa and Singh. dll HE] 


2. Calculations 

Introducing a different creation operator in each sublattice, A’*', H’*', C*” and 
D’t', the tight-binding term of the Hubbard Hamiltonian, is 
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where we set the hopping constant equal to I. 

We consider for now only the sublattice A (we add the other sublattice 
terms later on). The interaction term of the Hubbard Hamiltonian is, as usual, 
Hu = UYl We assume that the magnetic moments align in the 

r 

x-y plane, so that {Sz) = 0 and the Hartree term becomes where (n) 
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is the density of electrons on each sublattice (here assumed to be the same on 
all of them). 





(a) (n, U) phase diagram 


(b) (n, U) phase diagram with(c) (t/, fi) phase diagram with 
phase separation phase separation 



(d)m(n, t/) {e) qx{n,U) {i) qy{n,U) (g) E(n,U) (h) ii(n,U) 


Figure 2: (a) Mean-field phase diagram for the usual 2D Hubbard model: The system displays 
antiferromagnetism (AF), ferromagnetism (F), paramagnetism (P) or spiral phases {q / 0, tt). 
The antiferromagnetic state q = (7r,7r) only occurs for n = 1 (half-filling), (b) (n, U) and (c) 
{U^ii) mean-field phase diagrams for the 2D Hubbard model, allowing for phase separation 
(yellow region), (d) m, (e) qx-, (f) qy^ (g) Emf and (h) p, as functions of the doping n and 
Coulomb interaction t/, for the 100 x 100 2D Hubbard model. 


The Fock term includes averages like = {S^) = {Sax + whose 

values depend on the magnetic phase. Let us assume the average spin in the 
sublattice A is 


{Sta) = ^[cos(g-r:4),sin(g-r:4),0]. 


( 2 ) 


The vector q = {qx^Qy) defines the magnetic phase of the system. In /c-space we 
have 
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where Lu_c. is the number of unit cells, which gives 


while all the other mean values in the summation of Eq. vanish. The Fock 
term in Fourier space is 
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Adding the tight-binding, Hartree and Fock terms, the Hamiltonian Hmf reads, 
in the {Aj:, Bj:, C^, D^, basis, 


f Ht{k) Hm 
\ Hi Ht{k + 2^ 


plus the diagonal term 
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Here, Ht{k) is the tight-binding term (Eq. ^ of the Hamiltonian in fc-space. 
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is the diagonal matrix, = diag(Ay^, A^, A^, Aj:f) with 
UrriA 


Ad = —- 


Ac 


2 ’ 
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Ab = -^ 


Ad = 
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3. Results and discussion 


By setting tua = tub = me = mn = m^ we recover the MF magnetic phase 
diagram of the usual 2D Hubbard model, consistent with the ones obtained by 
several authors pTl 1131 ITS] for zero temperature, as presented in Fig. 2a In order 


to obtain such a diagram, one minimizes either the MF energy Emf using the 
electronic density n as an independent variable, or the thermodynamic potential 
^MF using the chemical potential /i, with respect to the site magnetization 
amplitude m and the order parameter q = {qx^Qy)- These parameters define the 
magnetic phase of the system. 

A solution with m = 0 is paramagnetic and is usually g-degenerate, while 
solutions for m 7^ 0 are in general unique. In the latter case, the wave vector q 
specifies the type of magnetic ordering. For instance, q = (0, 0) for the ferromag¬ 
netic phase, q= (tt, tt) for the antiferromagnetic phase and all other choices for 
spiral phases. In the example shown in Fig. we have qx = f/18 and qy = 7r/6. 
Additionally, in the same example, the magnetization amplitudes (denoted by 
the size of the arrows) are tua = mo = mi and mB = me = m 2 < mi. Com¬ 
paring, for each pair (n^U) or (t/,/r), the data obtained for m (Fig. [M[ ), qx 
(Fig. 2e) and qy (Fig. |^, the MF magnetic phase diagram displayed in Fig. 2a 
ensues. For some values of /i, there is more than one pair (g, m) which minimizes 
the thermodynamic potential. In those cases, a first-order phase transition in 
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(a) (n, U) phase diagram (b) (U, (i) phase diagram with phase 

separation 





(f) iJ.{n,U) 



Figure 3: Top figures: (a) Mean-field phase diagram for the 2D Hubbard model with sublat¬ 
tices with independent magnetization amplitudes. In the central (green) region of the phase 
diagram, configuration (i) is the one that minimizes the MF energy, while configuration (ii) is 
the most stable in the red regions, (b) (D, mean-field phase diagram with phase separation 
occurring on the borders of the green and red regions (thick solid lines). This phase diagram 
displays spiral and non-spiral ferrimagnetic phases (green regions) with fixed particle density, 
n = 0.5, refiecting the opening of an energy gap in the MF quasi-particle bands. Bottom 
figures: mi(n, f/) and m 2 {n,U) for the MF 2D Hubbard model with two sublattices in (c) 
configuration (i) and (d) configuration (ii). (e) Ground state MF energy for usual 2D (blue), 
case (i) (green) and case (ii) (red) for = 9 and U = 19. (f) as a function of the doping n 
and Coulomb interaction U. All figures represent numerical results obtained for the 100 x 100 
2D Hubbard model with two sublattices with different magnetization amplitudes. 


the order paranieters occurs. When using n as a basic variable (and posteriorly 
calculating /i = dEjdn ~ /S.Ej/\n using the data in Fig. 2g), n seems to be 
multiply defined for some values of which implies instability (e. g. of the 
spiral phase for U = 15). The use of /i as a basic variable solves this ambiguity 
and leads to plateaus in the chemical potential in the regions where 

phase separation (PS) occurs (see Fig. |2h[ ). In each PS region of the diagram, 
two spatially separated phases occur: the ones immediately to the left and to 
the right of the PS region in question (see Fig. 2b). The two phases have dif¬ 
ferent electronic densities, such that the electronic density of the whole system 
amounts to n. In Fig. we show the same phase diagram as in Fig. but 
using /i as the independent variable. The colors of corresponding regions are 
the same for easier reading. The thick solid line indicates a discontinuity in n. 
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In this work, the magnetic phase diagram for the Hubbard 2D model com¬ 
prising four sublattices is obtained by finding the magnetization amplitudes 
{mA^rriB^mc^mD) and the vector q which minimize the energy. We consider 
two situations: (i) ttia = ttijj = mi and ms = tuc = Tn 2 ] (h) mA = tuc = 


and mB = tud = ^ 2 - 

The ground state magnetization amplitude of the usual 2D Hubbard model 
is proportional to n for each value of U in the ferromagnetic phase (see Fig. 2d). 
When a spiral ferrimagnetic phase is allowed, it was found that, near zero filling 
(n = 0) and half-filling (n = 1), the ground state magnetization remains the 
same as in the usual 2D case (see Figs. 3c and 3d). This means that in these 


regions, the ground state magnetization is still constant throughout the whole 
lattice. However, as one moves to intermediate n, one finds that mi and m 2 
become distinct, as shown in Figs. 3c and 3d for cases (i) and (ii) respectively, 
where mi and m 2 are displayed as a function of n and U. These figures show 
two sheets reflecting the separation of the magnetization amplitudes. The colors 
green for case (i) and red for (ii) are used on all plots of Fig. For intermediate 
Ailing, the system is able to lower its energy by adopting different magnetization 
amplitudes on sublattices 1 and 2 in both cases (i) and (ii). This is shown in Fig. 
I^for U = 19 and U = 9. Depending on the region of the phase diagram one 
analyses, configuration (i) or (ii) may have the lowest energy, as shown in Fig. 
[3a| In this figure, we added another layer on top of the usual 2D MF magnetic 
phase diagram, showing which of the two-sublattice configurations considered 
has the lowest energy in the ferrimagnetic region: green for case (i) and red for 
case (ii). Furthermore, the energy was minimized with respect to Qx and 
while using the new magnetization values, but it was found that only very small 
changes in q occur, i.e., despite the changes in magnetization amplitudes, the 
magnetic phases in the diagram remain the same. For this reason, the magnetic 
phases are shown as being the same as those of the usual 2D model. 

The mean-field energy dispersion relation of the usual 2D Hubbard model 
displays two bands. Electrons occupy the lowest band until half-Ailing (n = 
1) and then proceed to occupying the higher band. As can be seen in Fig. 
|2h[ the fermionic density increases with the chemical potential until the phase 
separation region is reached. In this region, the chemical potential is constant 
despite any increase in the number of particles, up to half-Ailing. At this point, 
any increase in n induces a jump in the value of /i, equal to the energy separation 
between the two energy bands (called the energy gap). As the plot in Fig. 


2h 


only goes up to half-Alling, we see /i increasing smoothly until it reaches the 
phase separation region, followed by a plateau and a jump at n = 1. In both 
cases studied in this work, with the lattice divided into two sublattices, the 
energy bands open a gap at quarter Ailing (n = 0.5), as shown in Fig. 
Another gap appears at three quarter Ailing (n = 0.75), but only the phases 
with n < 1 are shown in Fig. |3b| The plot in Fig. is again the same as 
Fig. but using /i as the basic variable. In this diagram, the green region 
corresponds to n = 0.5, therefore only conAguration (i) for the spiral and non 
spiral ferrimagnetic phases is present in the phase diagram. The dashed line 
separates the ferrimagnetic region from the spiral ferrimagnetic one and the 
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thick solid lines denote again discontinuities in n. 


4. Conclusion 

Having addressed the possibility of a spiral ferrimagnetic phase in the mean- 
field phase diagram of the 2D Hubbard model, we conclude that, for interme¬ 
diate values of the interaction U and doping n, the spiral ferrimagnetic phase 
is the most stable phase in the (n, U) phase diagram. Higher values of U lead 
to non-spiral ferrimagnetic phases. We emphasize the case of intermediate n 
and higher [/, for which the ground state does not appear to be purely ferro¬ 
magnetic, contrasting with results by other authors. m Additionally, allowing 
phase separation and replacing n by /i as the independent variable, the (/i, U) 
phase diagram displays, in a considerable region, spiral (for intermediate values 
of U) and non-spiral (for higher values of U) ferrimagnetic phases with fixed 
particle density, n = 0.5, reflecting the opening of an energy gap in the mean- 
field quasi-particle bands. We further note that generalizing the ferrimagnetic 
phase to cases where more than two different magnetization amplitudes are al¬ 
lowed should lead to even more stable ferrimagnetic phases in certain regions 
of the phase diagram. Preliminary results with three different amplitudes are 
consistent with this conjecture. All these results provide strong evidence on 
the stability of the spiral ferrimagnetic phase in the mean field magnetic phase 
diagram of the 2D Hubbard model. 
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